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IN THREE DIMENSIONS WITH A SINGLE POINT SOURCE 



KLAUS FLEISCHMANN, CARL MUELLER, AND PASCAL VOGT 

Abstract. In a recent work, Fleischmann and Mueller (2004) showed the existence of a super-Brownian 
motion in R"*, d = 2,3, with extra birth at the origin. Their construction made use of an analytical 
approach based on the fundamental solution of the heat equation with a one point potential worked out 
by Albeverio et al. (1995). The present note addresses two properties of this measure-valued process in 
the three-dimensional case, namely the scaling of the process and the large scale behavior of its mean. 

1. Introduction 

A super-Brownian motion in R with a single point source 5q was constructed in Englander & Fleischmann 
|EFOO| . It was shown that its expected mass grows exponentially in time, and is in the mass-rescaled limit 
distributed in space as x i-^ e" -1^1. In Englander & Turaev |ET02j it is even proved that the random 
measures themselves grow in law exponentially as time increases, and are otherwise in the mass-rescaled 
limit spatially situated with the same shape except an overall random factor. The probabilistic effect 
behind the non-trivial existence of the model is the fact that a Brownian particle in M hits the origin with 
certainty and that it has there a non-degenerate local time, serving as an additional birth rate for the 
random creation of mass. 

In higher dimensions, a Brownian particle fails to hit the origin, and a local time would degenerate. 
Nevertheless, Fleischmann & Mueller FM04 succeeded in constructing a super-Brownian motion in R"^, 
d = 2, 3, with a single point source. They heavily used well-known analytical facts from mathematical 
physics concerning Laplace operators with one-point-potentials. Heuristically, some additional rescaling 
enters the regularization of the delta function (serving as single point source). Properties of this new 
super-Brownian motion are not known so far. The purpose of the present note is to get some progress by 
studying its scaling and the large scale behavior of its expectation in the three-dimensional case. 

1.1. The heat equation with one-point-potential. The Schrodinger equation with a one-point-poten- 
tial is studied in quantum theory to describe singular electromagnetic effects on quantum particles, see e.g. 
the monograph Albeverio et al. |AGHKH88I Part I] . By analytic continuation, solutions to the Schrodinger 
equation can be (at least formally) obtained via solutions of the heat equation. 
Formally, the heat equation with a one-point-potential is given by 

(1) = Au + 4"^u=: A^"^M, 

where dt denotes the derivative with respect to time, A is the d-dimensional Laplacian, and u : (0, oo) x 
M_i_ is a time-space field, where := M'* \ {0} with the Euclidean metric is locally compact. If we 
denote by B^{y) an open ball around y G M.'^ oi radius e > 0, then having in mind that e^'^lB^(o) ~ Sq , 
the operator A*^") := A -I- S^"^ is heuristically the limit as e | of the operator 

(2) 4") := A + /i(d,a,£)e-'^lB^(o), 

where h{d,a,s) is some additional rescaling factor which depends on a parameter a at least. Restricting 
to d = 3, the function h can be chosen as 

(3) /i(3,a,£) := ^e-STT^ae^ a G M, e > 0, 
(cf. IAGHKH8SI (H.74)]). 

Physically, a in the case a < is related to the scattering length sIq '.— — (47ra)~^ of the free Laplace 
operator A with respect to the interaction Laplacian A*^"'. Roughly speaking, the scattering length de- 
scribes the average distance a free particle manages to go before any interaction takes place. So, if a J, — oo 
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the scattering length sIq. J, becomes smaller and we expect more interaction. For a > there is no proper 
physical interpretation of sIq, as the point spectrum of A^"^ is empty (see (AGHKH88i Theorem 1.1.4]). 
The fundamental solution to the equation 

(4) atM = A("'u on (0, oo) X M'', d = 2,3, 

which provides the basis for the analytical construction of the superprocess in |FM04| . have been computed 
in Albeverio et al. [ ABD95j . In d = 3 (the two-dimensional case is analytically more delicate, which is the 
reason we restrict to d = 3), the one-point-interaction heat kernel p" is given by 

(5) p'^{x,y)=pt{x,y) + --—pt{\x\ + \y\)~--— e-*"""pt(u + |a;| + jyl), 

< > 0, x,y ^ 0, where p is the usual free heat kernel defined by, 

(6) Pt(x, y) := {Ant)-"^^ exp ( - |y - xl^At), 
and with a slight abuse of notation, 

(7) pt(r) (47rt)-'^/2exp(-rV4i), t>0, r>0. 
Also recall the scaling of the free heat kernel, i.e. for all k,t > and x,y £ M'^, 

(8) ptix,y) = k^^'pkt{k'/^x,k'^'y). 

Note, that the last term in © is always finite and disappears for a = 0. Moreover, a i— > is pointwise 
continuous and decreasing, and we have the (pointwise) convergences t +oo as a J, — oo (i.e. the 
fundamental solution explodes which can be interpreted as immediate interaction), whereas p" j p as 
a t +CXD leads the free case (i.e. the interaction disappears). 

Rigorously, the family {A^"' : a € K} of operators are defined as all self-adjoint extensions on the 
Hilbert space £.'^{W\ dx) of the Laplacian A acting on C^^iK^), the space of unboundedly differentiable 
functions on R'' = R'' \ {0} with compact support (see e.g. |AGHKH88I Chapters I.l and 1.5]). Hence, 
although the p" differ from the free heat kernel p, they solve the heat equation 

(9) a(P«(x,y) = Ap«(x,y) on (0, ^) x R^ 

with the Laplacian A acting either on the variable x or y. In particular, {t,x,y) »— > pf{x,y) is jointly 
continuous on (0, oo) x R'^ x R'^. Let us denote by S" the semigroup associated with the kernel p", i.e. 

(10) S^^ix) / dy^iy)pfix,y). 

Note that 5" is not a contraction semigroup and so there is no stochastic process generated by this flow. 
The following Lemma shows that the kernel p" has a similar scaling behavior as the free heat kernel p. 

Lemma 1 (Scaling of the p"). We have, for all k,t > 0, x,y E M^, and a £ M, 

(11) pnx,y)^k^^'pi;''''^{k'f'x,k'^'y). 

Proof. That follows immediately from the definition 10) of the p" and the scaling of the free heat kernel 
p. □ 

1.2. The flow associated with the one-point-interaction heat kernel. This section is devoted to 
introduce a space of functions $ on which the flow S" acts as a strongly continuous linear semigroup (see 
|FM04[ Section 2] for details). Let </> denote the weight and reference function 

(12) (j){x) -.^ \x\-\ x e = \ {0}. 

For fixed g G (1, 2), let Ti — TL^ denote the space of measurable functions on R'^ for which 

(13) y\\^:=ff dx^{x)\^{x)\<^y^' <cx. 

Then (7i, || • is a Banach space, where as usual we do not distinguish between equivalence classes and 
their representatives. Now, let $ — denote the set of all continuous functions : R'^ ^ M such that 
ip £ H and 

(14) < (y5 < C for some constant C — C^p > 0. 

We endow $ with the topology inherited from Ti. Note that the set C+„„ = C(^„j(R'^) of all non-negative, 
continuous functions on R'^ with compact support is contained in $. We remark that (p £ ^S? might have a 
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singularity at a; = of order \x\ ^ with < ^ < 1 . The Hnear semigroup 3°^ introduced in (|10|l is strongly 
continuous on the cone $ = <I>^, cf. Corollary 2.12 in FM04 . 



(16) 



1.3. Super-Brownian motion with a single point source. Denote by = M.{M.'^) the set of all 

(Radon) measures /i on R'^ such that (/i, (^s) < oo for all (ys G $. Recalling that C+,i„ C endow M with 
the vague topology. 

Fix a constant 77 > (branching rate) . Suppose < /3 < 1 (the finite variance branching case (3=1 has 
been excluded in |FM04| for d = 3 for technical reasons) . Then for each a G R, there is a non-degenerate 
A^-valued (time-homogeneous) Markov process X" such that for (deterministic) starting measures /i G 
Ai and for G $, 

(15) -logE^exp(X,",-(^) = (^,t-(t,.)), t>Q, 

where {w(i, x) : t > 0, a; G R'^} is the unique non-negative solution of the integral equation related to the 
$- valued evolution equation 

J = A(")w - ?7w^+'' on(0, cx)), 
(see jFM04l Theorem 4.4]). That is, 

(17) v{t,x) = / dy p'^{x,y)(p{y) -ri ds / dy p'^_,{x,y) v'^+'^ {s,y), 

Jt^ Jo Jii^ 

t > 0, X G R'^. Clearly, the first moments of X" are determined by the S*" flow to be 

(18) E^(xr,<^) = (/.,5,», 

for all starting measures fi G M, t > 0, and G 

2. Large scale behavior 

2.1. A scaling property. 

Proposition 2 (A scaling property). Let t,k> 0, /i G A4, and a,Xk G M. Then 

(19) {fc-i/^X,^,-"(fci/2.)|Xo^'=" = fc^/''A^(fc-i/2.)} = {xf'''^*" I Xo^'''^-" = . 

Of course, the cases At = k^^^^ or even a = are particularly nice. 
Proof. For G $ fixed, 

(20) {k-^/^X^rik'/'dy),^) = {X^r^k^'^'^ik-'^'-)), 
hence, by and ((T7jl . 

-logEfci/,^(,-V2.)exp(fc-i//'x,^^'="(fci/2rfy), ^) = -logE,:/^^(fc-:/2.)exp(X,^r, fc^i/Mfc^'^' ' )> 

(21) = (fci/^^(fc-i/2.)^^(fc^^ ^ {n,k^/^vikt,k^/^-)) 

where {w(i',a;') : t' > 0, x' G R'^} is the non-negative solution of the integral equation related to the 
function-valued evolution equation 

( dtv ^ A^^''°'^v -r]v'^+'^ on(0,oo), 

^^^^ \ w(o+, ■) = fc"'/Mfc~'/'-)- 

More precisely, 

(23) k^/l^v{kt,k^''^x) ^k^/l^ j dyplf{k^''^x,y)k-^'Pip{k'^''^y) 

- k'/^v /'*ds / dyptt^^ik'/'x,y)v'+Pis,y). 
Jo Jm 

By a change of variable, 

(24) k^^f^v{kt,k^^^x) = / dy fc3/2p^^^"(fci/2:E,A:i/2y)^(y) 



/ dsk dye/^pll^^^{k^/^x,k^'^y)v^+f^{ks,k^'^y). 
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Hence, by Lemma ^ 

k^/l^vikt.k^'^x) = I dy pf^'^'''^{x,y)ip{y) - k^^/^rj f ds k f dy pf_!l^''°'{x,y)v'^+P{ks,k'^'^y). 
Since 1 — (1//^)(1 + /3) = wc sec that k^/Pv{kt, k^^'^x) —: Wk{t,x) satisfies the equation 

(25) Wk{t',x') = / dy p;,' ^""'{x ,y)ip{y) -ri / ds / ^^^'^(s, y), 



t' > x' <E W' . By uniqueness of soiutions of the log-Laplace equation (|17|l and by \ib\ . claim p9|) 
follows. □ 

2.2. Expectation of the scaled X". Before we can state the result, we have to introduce some notation. 
The limiting measure will be expressed by means of the kernel 

(26) '^"(^'^^^-RM^'^'^'^-Rm/o '^--"">*("+M)' 

for a e R, t > 0, and x,y € R'^. Note that the integral is always finite, hence for a — the second term 
disappears. Moreover, the kernel is always non-negative. This holds trivially whenever a < 0, and to 
see this for a > 0, use the estimate 

(27) Pt{u + \y\) <pt{\y\). 
We extend the definition of by setting 



(28) ^?{x,y)-= 



_ j 0, if a — +0O, 



oo, if a = — oo. 



The so defined kernels i?" turn out to be pointwise continuous in a G [— oo, +oo] (which follows from the 
arguments of the proof of Theorem |31 below). 

Theorem 3 (Large scale behavior of the mean). Fort > 0, a,Xk £ R, and all starting measures 
Xq = /i G satisfying (/i, 0) < oo, we have the convergence in Ai, 

(29) hmfc-i/2E4x,^r(fc'^'rf2/)] = {l^,'df{-,y))dy, 

k\oo 

provided that a* := limfc|oo k^^^XkCt G [~oo, -l-oo]. 

Proof. Fix (f G C+„„(R'^). Using formula (|18() for the first moment of and substitution, we obtain 

k-'/'E,{X^r,^ik-^^'-)) = k-'^' f l^{dx) ( dyplri^,y)^{k-'/'y) 

(30) =k-'/^ f f^idx) f dyk^'^plr{^,k^'^y)v{y)- 



By Lemma n] this is equal to 

(31) / ^{dx) I dypf'^-^{k-''^x,y)^{y). 



Inserting according to definition |5j of p", we get three terms, we will deal with separately. 
1° [First term). The first term equals, 

(32) / ^Ji{dx) / dypt{k-^'^x,y)^{y). 



This double integral is finite and vanishes as fc t oo. To see this, let us restrict the outer integral first to 
|a:| > K where we specify K >1 later. We call this restricted integral Ik ■ We use Lp < Ccj) (since i^s G $) 
and, with S denoting the free heat flow, 

(33) St(t) <C(I), t> 0, 
with changed constant C (see |FM04I Lemma 2.1]) to arrive at 

(34) Ik < Ck-^/^ [ fi{dx) (j){k-^/^x) ^ C [ n{dx)(j){x) 

J\k-^/^x\>K J\x\>k^/'^K 
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which can be made arbitrarily small uniformly in k by choosing K sufficiently large (by our assumption on 
n). It remains to deal with the case |a;| < K for fixed K. We split the internal integral in ()32|l as follows. 

First, if |fc-i/2^ -y\> \y\/2, then 

(35) Mk-'/^x,y)<pt{\y\/2), 
which leads to the bound 

(36) fc~i/2 /■ f,{dx) f dypt{\y\/2)ipiy) ^0 as fc t oo, 

J\x\<K JR3 

the /i((ix)-integral is finite as (/i, cp) < oo. On the other hand, if \k~^^^x — y\ < \y\/2, then — fc^-'^''^|a;| + < 
j/l/2 which implies |y| < 2fc~^/^|a;| < 2K. Hence as pt{k~^^^x,y) < Ct^^^"^ and (p < Ccj), we get the upper 
estimate 

(37) Ctk-^/^ [ iJ-{dx) I dy(f>{y) — >0 as fc T oo 

J\x\<K J\y\<2K 

(the dy-integral is finite, since we are in dimension three). 
2° (Second term). The second term reads 

(38) / fi{dx) f dy -^pt{k-'/'\x\ + \y\)^{y) Ilk. 

Jii^ Jr3 \x\ \y\ 

Observe that, 

(39) Pt{k~^^^\x\ + \y\) ^ pt{\y\) as fc T oo. 
We can apply the monotone convergence theorem to obtain the limit, 

f 2t f 1 

(40) jim Ilk = / nidx)— dy —pt{y)<fi{y), 

^^Too \x\ 7r3 \y\ 

where finiteness follows from (/i, </!>) < oo and Lp < Ccj). Hence, this summand gives the first part of the 
kernel i?" . 

3° ( Third term) . It remains to insert the scaled third term from (|SJ| into (|31ll which reads as 



(41) fc-V^ / ^i{dx) / dy ° ,7" / due-^^^' ^''^^p^{u + k-''M + \y\)v{v)■ 
We distinguish several cases: If Afea — for all sufficiently large fc, then the third term disappears and we 
are done. From now on assume A^o: ^ for all fc. Substituting u i— > 47rk^/'^\\ka\u into (|41|l yields, 

(42) / Mrf-)/ ^^ -2^Bign(A.a) T ,-.ign(p..) / J fc-i/2|^| ^ (^)^ 
Jr3 Jk3 \x\\y\ Jq V47rfci/2|Aj,a| / 

Now let k^^^\Xk a\ — > 00. We may consider a monotone subsequence of k^^^XkCt. Clearly, 

(43) p^i^ ^^^J^^^^^^ +k-^/2\^\ + \y\)^p,{\y\) as fc T 00, 

and by monotone convergence the expression (|42|l converges along the subsequence to 

(44) -/ ^,{dx)^^^^^^ r due-^'s^^»'> f dy ^p,iy)ipiy), 



Fi Jo \y\ 

which is independent of the choice of the subsequence. Note that 

1 if sign(a*) = 1, 
hoo if sign(a*) = — 1. 



/"OO 

(45) sign(a*) / e-^'s"^"') = 

Jo 



/2|Afea| 



In the first case the second and the third limiting terms cancel. 
Next, let fci/^AfeO; 0. Note, that 

In this case the double integral (|42|) is in absolute value bounded by 

(47) / f^idx) I dyj^ e-\y\''''^{y) r due-p,( J ), 

7r3 \x\\y\ Jo V47rfci/2|Ai,a|/ 
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which tends to as fc f cx) as the fi{dx) and dy-integrals are finite and in the dw-integral the pt-term 
compensates the e". 

It remains to deal with the case /c^/^A^a a* £ W- . Note, that we only have to justify to change 
the limit and integration in (|42|l . as substituting u ^—^ (47r|Q;*|) u leads to the desired expression. To 
justify the interchange, we estimate as in (|46|) . The resulting fi{dx) and dy-integrals are independent of k 
and finite, whereas to dominate in the second integral we use fc^/^|Afca| < |a*| + 1 for all sufficiently large 
k. □ 

Remark 4 (Large scale total mass). Taking Xk = I and choosing formally = 1 as test function in 
yields 

if a > 0, 

(48) limk-'/^E^{X^„l) = <! 2t{fi,cj)) [ dy ^ pt{y) if a = 0, 




feToo I Jr3 \y\ 

if a < 0. 

A rigorous argument can be given along the lines of the previous proof. O 

2.3. Discussion and open problems. Let us comment on the three cases a* — +oo, a* G R, and 
a* = — oo in Theorem|3 In the first case, the limiting mass disappears, more precisely, the scaled expression 
E^[X^/="(fci/2 jy)] is of order o(fci/2). Roughly speaking, if a* = +0O, then there are no interactions in the 
scaling limit (free case). In the second case, a* G M, the former expectation is about fc^/^ (/i, ( • , y)) dy. 
Note that these measures are decreasing in a* . Finally, if a* — ~oo, we have immediate interaction in the 
large scale limit leading to the explosion of the expected mass. 

Clearly, to describe only the large scale behavior of the expected processes is unsatisfactory. It is desirable 
to get something similar for the processes themselves. Recall that in the one-dimensional case the large 
time behavior of the process itself is known from (ET02| . However, we stress the fact, that the process in 
three dimensions is expected to have quite different features. For instance, if a = 0, then according to 
Remark 0] the total mass grows with a power order, whereas in one dimension the growth is exponential. 
Moreover, in the three-dimensional case one needs additionally to contract the normalized measures to get 
a limit. For the measures themselves, scaled as in Theorem |21 there might be extinction in law despite 
convergence of their expectations. 

Another open problem is the large scale behavior of EX" in the two-dimensional case, in which the 
fundamental solutions p" from |ABD95| are analytically more delicate, see e.g. (FM04i formula (2.30)]. In 
particular, a scaling property as in Lemma ^ is not available. 
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